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The method of asymptotic expansions is used to build an approximation scheme relevant to 
celestial mechanics in relativistic theories of gravitation. A scalar theory is considered, both as a 
simple example and for its own sake. This theory is summarized, then the relevant boundary 
problem is seen to be the full initial-value problem. It is shown that, with any given system of 
gravitating bodies, one may associate a one-parameter family of similar systems, the parameter 
measuring the gravitational field- strength. After a specific change of units, the derivation of 
asymptotic expansions becomes straightforward. Two hypotheses could be made as to which 
time variable has to be used in the expansion. The first one leads to an "asymptotic" post- 
Newtonian approximation (PNA) with instantaneous propagation, differing from the standard 
PNA in that, in the asymptotic PNA, all fields are expanded. The second hypothese could lead to 
an "asymptotic" post-Minkowskian approximation (PMA) allowing to describe propagation 
effects, but it is not compatible with the Newtonian limit. It is shown that the standard PNA is 
not compatible with the application of the usual method of asymptotic expansions as envisaged 
here. 



1. Introduction 

In order to evaluate the predictions of relativistic theories of gravitation, e.g. for the motion of 
celestial bodies in the solar system, it is necessary to use some approximation method. In 
Einstein's general relativity (GR), the most widely used method for that purpose is the so-called 
"post-Newtonian approximation" (PNA), elaborated notably in Refs. [1-5]. The PNA is the 
method used to consolidate one of the key tests of GR, namely that it predicts so accurately 
Mercury's residual advance in perihelion. The textbook derivation of that prediction is based on 
the ideal situation in which Mercury would be a test particle orbiting in the static, spherically 
symmetric gravitational field of the isolated Sun. Using the PNA, one finds that the same 
residual advance is predicted if one assumes a more realistic situation where, in particular, the 
perturbations due to the other planets are accounted for [1, 3, 5]. More generally, the aim of the 
PNA was to build a consistent relativistic celestial mechanics and to compare its predictions with 
astronomical observations. 

Basically, the standard PNA consists in defining a small parameter £ that characterizes the 
weakness of the gravitational field [1, 4], in expanding the gravitational field in powers of £, and 
in deriving approximate equations from this expansion and the exact field equations. A first 
question arises: in which sense does the solution of the approximate equations actually approach 
the relevant solution of the exact equations? This question is complicated, in GR, by the special 
structure of that theory, in particular by the necessity of imposing a gauge condition [6-7]. To the 
author's knowledge, this question has not received a clear answer today. A rigorous derivation of 
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a PNA scheme, resulting from the assumed existence of a one-parameter family of metrics and 
mass tensors obeying the Einstein equations, has yet been proposed by Rendall [8]. But, as noted 
by Rendall, his scheme differs from the standard PNA in that, in his scheme, all variables, thus 
gravitational fields as well as matter fields, are expanded. 1 In contrast, in the standard PNA, only 
the gravitational field is expanded (the leading expansion being that of the space-time metric y, 
which gives the expansions of the Christoffel symbols and the Ricci tensor). As a consequence, 
Rendall' s work [8] does not seem to provide a direct justification for the standard PNA. It rather 
gives motivation for developing the application of the method of ordinary asymptotic expansion, 
even if it leads to difficulties in the comparison with observations [12]. Moreover, Rendall [8] 
does not attempt to describe how the one-parameter family of metrics, that he considers, could be 
related to a family of physically relevant boundary value problems. 

A second difficulty is that the approximate equations for the gravitational field, obtained by 
using the PNA, turn out to be Poisson equations for the Newtonian potential and some 
additional, "post-Newtonian" (PN) potentials. Thus, whereas an essential new feature of 
relativistic theories of gravitation is that gravitation propagates with speed c, it is found that 
gravitation propagates instantaneously at the PNA; this is not only true at the first (1/c ) 
iteration 3 , but also at the second (1/c 4 ) one [13]: only at the "2Vi" (1/c 5 ) iteration of the PNA does 
one find out propagation effects, including reaction effects to the emission of gravitational waves 
[14]. But, as noted by Ehlers et al. [6], the determination of the lowest-order time-odd terms in 
the expansion of the space-time metric, which is crucial in the derivation of the 2Vi iteration of 
the PNA, was not achieved in a self-consistent way in Ref. 14: this determination used results 
obtained with a different approximation method, namely the standard linearized theory of 
gravitational radiation; in brief, the statement of the IVi PNA in Ref. 14 lacks the straightforward 
character that one might a priori expect from an iterative approximation scheme. The "post- 
Minkowskian" approximation (PMA) schemes differ from the PNA in that the approximate 
equations for the gravitational field in the PMA involve a d'Alembert (wave) operator in the 
place of Laplace's operator, thus automatically maintaining the effects of propagation [15-16]. 
Damour & Schmidt [17] gave existence proofs of one-parameter families of metrics solutions of 
the Einstein equations, and showed that the successive coefficients in the Taylor expansion (with 
respect to the parameter) of a such family satisfy the hierarchy of equations encountered in the 
successive iterations of the PMA. As for Ref. 8 (devoted to the PNA), they do not link the 
appearance of their family of metrics with a family of physically relevant boundary value 
problems. Moreover, they study Einstein's vacuum field equations, whereas celestial mechanics 
should consider a system of extended bodies [9]. In vacuo, there is only the gravitational field to 
be expanded, hence the work [17] does not apply directly when the standard PMA, consisting in 
expanding the gravitational field only [9, 15-16], is used [9] for extended bodies. 

The aim of the present paper is to build a general approximation scheme based on the 
usual method of asymptotic expansions and to obtain the explicit expansions of the fields, 
equations and boundary conditions up to the order 1/c 2 included. Preliminary results have been 



1 One must expand all unknown fields if one searches for asymptotic expansions, in the usual sense, of solutions to 
systems of partial differential equations. Typical examples can be found in transonic gas dynamics [10] and in the 

homogenization theory as applied, e.g., to porous media or to composite materials [11]. 

2 

A kind of expansion of some matter variables such as the components of the energy-momentum tensor T is used in 
the standard PNA, but this is a "composite expansion" ; it occurs as a by-product of the expansion of y (due to the 

fact that any particular expression of T, e.g. that for a perfect fluid, does involve the metric y). 

3 

In appropriate units, one may assume that the small parameter e is simply 1/c with c the velocity of light, as is 
often formally assumed [see § 5.2]. 



2 



presented in Ref. 18. In this whole work, we shall consider a scalar, preferred-frame theory of 
gravitation, proposed in Refs. 19-20. (Cf. Ref. 21 for a review including recent results.) This 
choice has the following advantages: (i) in the scalar theory, no gauge condition is needed, so 
that the comparison between the proposed "asymptotic" scheme and the standard PNA will be 
unbiased; (ii) since the scalar theory assumes a "prior-geometrical" flat metric and since its exact 
field equation may be formulated in terms of this flat metric, the "flat space-time operators" 
needed to obtain a tractable approximation scheme occur naturally and do not raise any difficulty 
of interpretation, in contrast to the situation in GR [6-7]; (iii) the scalar theory is mathematically 
much simpler than GR, in particular the derivation of asymptotic expansions is very easy. Thus, 
the choice of this theory provides a simple illustration of the application of the usual method of 
asymptotic expansion and of its difference with the standard PNA and PMA. Moreover, this 
theory predicts just the same effects of gravitation on light rays as does GR, thus in the 
gravitational field of a static spherical body, and this is also true if, as one expects, the body is in 
a translation with respect to the preferred frame of the theory [22] . Due to the importance of this 
first experimental check, it is worth to test the scalar theory further and, for that goal, to develop 
a celestial mechanics according to that theory; the present paper may also be considered as a step 
in this direction. 

The plan is as follows: in Section 2, we summarize the scalar theory which is considered. In 
Section 3, we derive the exact equations for a perfect fluid. In Section 4, we show that the 
relevant boundary-value problem is the full initial-value problem. In Section 5, we describe how 
to associate, with the given gravitating system S, a one-parameter family (S e ) of similar systems. 
In our opinion, this is needed in order to define a meaningful asymptotic approach in a physically 
relevant situation: if we do not have such family at our disposal, it is rather in a formal sense that 
we may use an asymptotic method for solving the given physical problem. Of course, one is often 
led to use formal asymptotics in physics, but the presentation of a framework really suitable for 
an asymptotic method seems to be a progress. In Section 6, we derive the asymptotic expansion 
of the field equation and the local equations of motion for a perfect fluid, using a PNA and 
expanding all field variables, and we compare the equations thus obtained, with the equations 
obtained when one expands the gravitational field only. In Section 7, we discuss the boundary 
conditions for the expanded fields. We finish the paper by the Conclusion (Section 8) and a 
bibliographical note (Section 9). 



2. A short summary of the scalar theory 

Space-time is assumed to be equipped with both a flat metric 7° and a curved metric y, the latter 
being more directly related to physical measurements (the "effective metric"). The only equations 
of this Section that will be used in the remainder are those defining the effective metric in terms 
of the scalar field / [Eqs.(2.3) and (2.5)], the field equation (2.6) for /, and the dynamical 
equations (2.9)-(2.10). Obviously, they are independent and may be taken axiomatically, 
independently of the physical status that one wishes to attribute to that theory. The preferred 
frame E is assumed to be an inertial frame for the flat metric. By this, we mean that there are 
Galilean coordinates (xf) for y° [i.e., in that coordinates, (^v) = (Vnv) = diag(l, -1, -1, -1)], that 
are adapted to the frame E. ("Adapted coordinates" are such that any particle bound to the given 
frame has constant space coordinates [23].) The inertial time T = x°/c in the preferred frame is 
called the "absolute time". The scalar field of the theory is the contraction factor jS affecting the 
time interval dt x measured by a clock fixed at point x bound to the frame E (due to the dilation of 
the clock period in the gravitational field), or the square / = fi 2 . The latter defines the /oo 
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component of the physical space-time metric y in any coordinates (y), though with y Q = cT, that 
are adapted to the frame E: 



d; x /dr= J(y 00 ) F =4f= p. (2.1) 



The field/is also a potential for the assumed gravity acceleration g [20]: 

2 
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g = -yV/, (v/)' s y°f/j, (2-2) 

where (g° ij ) is the inverse matrix of the component matrix (g°i/) of g°, with g° the spatial part of 
the flat metric y° in the frame E (thus g° is an invariable Euclidean metric defined on the 
"preferred reference body" M, i.e. the set of the particles bound to E). The field jS still defines the 
dilation of an infinitesimal distance dl measured in the direction of vector g, as compared with 
the distance d/° evaluated using the Euclidean metric: dl = dl°/(5, while the infinitesimal distances 
in directions perpendicular to g are unaffected. (This makes the physical space metric g in the 
frame E discontinuous at the isolated points where g = 0. But g is merely the effective metric 
and, moreover, it remains bounded in the neighborhood of any such point: no theoretical 
difficulty occurs [22].) This gives the following expression for metric g in the frame E (in any 
coordinates (y /i ) adapted to E) : 



9ij = 9°ij + 



1 

1/ j 



h,, h n = ■'■ iJj (2.3) 
9 f, k fj 



The space-time metric writes, in any coordinates (y M ) adapted to E and such that y = cT : 

ds 2 =y^dy^dy v =f{d y y-dl\ dl 2 = g ij dy i dy i , (2.4) 

which gives 

7bo=f, Ya = -9ip 7o/ = 0. (2.5) 

Thus, the 70/ components are always zero in the ether frame E. This is the expression of the 
assumption that this frame admits a global simultaneity defined with the absolute time T. Except 
for cosmological problems, the field equation [19] may be written as [20]: 



1 

A/- — 
/ 



ffi 

u 



871 G n 

—t-G (y Q = cT), (2.6) 

c 



with A the usual Laplace operator, defined with the Euclidean metric g , and where G is 
Newton's gravitation constant, and 0= {T m )^ is the mass-energy density in the ether frame. This 
equation, which is merely space-covariant, implies that gravitation propagates with the velocity 
of light c (as measured with physical standards affected by gravitation). 

Motion is governed by an extension of Newton's second law [20]: 



4 Latin indices run from 1 to 3 (spatial indices), Greek indices from to 3. 
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F + m(v)g = DP/Dt x . 



(2.7) 



With the gravity acceleration g assumed (in the preferred frame E) in the theory [Eq.(2.2)], it 
implies Einstein's motion along geodesies of metric y, but only for a static gravitational field; 
however, Newton's second law (2.7) may be defined in any reference frame, and also for metric 
theories: for those, one adds to the g of Eq. (2.2) a certain velocity-dependent term, which gives 
geodesic motion in the general case [24]. In Eq.(2.7), Fo is the non- gravitational force, v is the 
modulus of the velocity v of the test particle (relative to the considered arbitrary frame F). The 
velocity v is measured with the local time t x synchronized along the given trajectory [25] (t x is 
given by Eq.(2.1) if the preferred frame of the present theory is considered, thus if F = E) and its 
modulus v is defined with the space metric h in the frame F (thus h = g given by Eq.(2.3) if F = 
E): 

v' = dv ! /d; x , VE[h(y,y)] 1/2 = (^v'v J ) 1/2 ; (2.8) 

m(v) = m(0).y v = m(0).(l - v 2 /c 2 )~ 1/2 is the relativistic inertial mass, P = m(v) v is the momentum, 
and DIDt x is the derivative of a spatial vector appropriate to the case where the Riemannian 
spatial metric h varies with time. In particular, with this derivative, Leibniz' rule for the time 
derivative of a scalar product v.w = h(v,w) is satisfied [20, 24]. 

Equations (2.2) and (2.7) define Newton's second law in the preferred frame, for any mass 
particle. It is thereby defined also for a dust, since dust is a continuum made of coherently 
moving, non-interacting particles, each of which conserves its rest mass. It then mathematically 
implies, independently of the assumed form for the space-time metric y provided it satisfies )fo = 
in the preferred frame, the following dynamical equation for the dust, in terms of its mass 
tensor T^ v = p*U^U v (with p* the proper rest- mass density and C/ 1 " = dy^/ds the 4-velocity) [26]: 

T^. v = V 1 - (2.9) 

Here b M is defined by 

b^) = \g jk fiT ]k , Ws-^oT- *. (2.10) 

[Indices are raised and lowered with metric y, unless mentioned otherwise. Semicolon means 
covariant differentiation using the Levi-Civita connection associated with metric y. The 
derivation of Eqs.(2.9)-(2.10) [26] is outlined in Ref. 21; they are valid in any coordinates (j/ 1 ) 
that are adapted to the frame E and such that y° = 0(7) with T the absolute time.] Equation (2.9), 
with the definition (2.10), is assumed to hold for any material continuum: accounting for the 
mass-energy equivalence, this is the expression of the universality of gravitation. 

3. Exact dynamical equations for a perfect fluid 

With the help of the identity 




(3.1) 



[in which the T ' fa 's are the Christoffel symbols of metric y], selecting coordinates such that y= 
det(7 MV ) = -1, and using Eq.(2.5), one may generally rewrite the spatial part of Eq.(2.9) as 
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r v , v + V jk T jk + \g lJ {g jkfi T Qk + fj T 00 ) = 0, (3.2) 

in which the T l j k 's are the Christoffel symbols of metric g and (gP) = (flfy) - . Note that 7= -1 in 

Cartesian coordinates for the Euclidean metric g° {i.e., coordinates in which g ^ = 5y), since 
Eq.(2.5) implies that J = fg with g = det(g,y) and since the assumed relation between g° and the 
physical space metric g implies that 

g = g°/f (flP = det(g<V)). (3.3) 

Let us apply Eq.(3.2) to a perfect fluid. The mass tensor is then 

T^=(p*+p/c 2 ) U^U v -(p/c 2 ) J^, (3.4) 

where /j* is the volume density of the rest- mass plus internal energy in the proper frame, 
expressed in mass units, 



We note that 



jl* = p*{\ +n/c 2 ). (3.5) 



^k9 jk = 9 il 9ikj9 jk -\9 U 9jk,i9 jk (3-6) 

and that, due to Eq.(3.3), one has 

9jk,i9 - g - g0 f - (3-7) 

Using this, it is straightforward to check that, for a perfect fluid, Eq.(3.2) takes the form 

d T {\iru i ) + d j {yu i u i ) + T) k yu i u k + t i k \iru k = y/fg i -g ij p J (3.8) 
(valid in Cartesian coordinates). Here 

fk ^\^d T g jk , ¥ =o+p/(c 2 f). (3.9) 

As to the "time" part of Eq.(2.9), one gets by lowering index /j in (2.9) and by writing the ji 
= component in Cartesian coordinates: 

d T (xi/f) + dj(iiffu j ) = (G/2)d T f + (l/c 2 )d T p. (3.10) 
4. The question of the relevant boundary value problem for a perfect fluid 
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To be definite, we assume henceforth that any matter is in the form of barotropic perfect fluids 
(one fluid per astronomical body), which means that /J,* in Eq.(3.5), hence also p* and 77, depend 
only on the pressure p, with 77 being given by [1] 



nip) = 1; 



dq P (4.1) 



p*(q) p*(pY 

(The integral must be finite.) Note that the source of the gravitational field [Eq.(2.6)] is given by 



o=(T ) E = 



P A 



V c J 



Yr P (4.2) 



and is hence determined by p, f and u. The field equation (2.6), together with the dynamical 
equations (3.8) and (3.10), make a system of five independent partial differential equations for 
the five independent unknowns p,/and u (or y/,/and u; or (7,/and u). All of these equations are 
nonlinear. The dynamical equations (3.8) and (3.10) make a nonlinear first-order system of the 
same nature as, though more complicated than, the nonlinear system constituted by the continuity 
equation plus Euler's equation for a perfect fluid in Newtonian gravity (NG) [see Eqs.(5.2) and 
(5.3) below]. Already in NG, it is not possible to consider that the gravitational field (the 
potential U^, in NG) is given, i.e. U^(x,T) known for all x and T ; for, in that case, Poisson's 
equation first determines the density p, and then Euler's equation alone [though with an initial 
data u(x,r=0)] determines u, so that the continuity equation has little chance to be satisfied. A 
similar situation occurs in the new scalar theory, as one easily checks. 

However, if the matter source cr(x, 7) is considered given, Eq.(2.6) is a quasi-linear second- 
order equation for the scalar gravitational field/. It may be written, in Cartesian coordinates, as 

\vf,n,v + B = 0, Aoo = l/f\ A Qi = At o = 0, A i} = -5 ih B = SnGa/c 2 -f^/f. (4.3) 

It is thus a hyperbolic equation, the characteristic equation to be satisfied on a characteristic 
hypersurface 0(jc ^ ) = being [27] A^ <j>^ <j) >v =0, i. e. here 



2 rl 

c f 



= (V0) 2 . (4.4) 



The hypersurface T = (or T = Const), i.e. 0Oc^) = x - Const = 0, cannot be characteristic, for 
this would mean that l//(x,0) = 0. Therefore, the Cauchy problem, in which one assumes an 
arbitrary initial data/(x,T=0) and djf(x,T=0) for all x, is very plausibly a well-posed problem, and 
this should be the case for a generic source a= o(x,7). On the other hand, in NG, one also may 
consider the source p(x, 7) as given: but Poisson's equation is an elliptic linear second-order 
equation for U^, so the problem is well-posed when one assumes, at any time T, a spatial 
boundary condition for ; typically the condition at infinity, C/n = 0(l/r) and V£/n = 0(l/r ) as 
r — > oo, is assumed. 

Now, in NG, the natural boundary value problem for the full system of equations is a 
mixed initial/ spatial-boundary problem, in which one assumes a spatial boundary condition at 
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any time for U, as just mentioned, plus an initial data for the velocity u and the density p (or the 
pressure p). Since (i) the dynamical equations (3.8) and (3.10) have the same structure as the 
system made of the continuity equation plus Euler's equation in NG, and (ii) the natural 
boundary condition for/ is an initial condition, whereas the natural boundary condition for C/n is 
a spatial boundary condition (due to the hyperbolic nature of Eq. (2.6), as opposed to the elliptic 
nature of Poisson's equation), we are led to consider the full initial-value problem as the relevant 
boundary value problem in the scalar theory. That is, one should assume the data cr(x,0) [or 
p(x,0), or p*(x,0)], u(x,0),/(x,0) and djf(x,0), for all x. 

5. The situation suitable to derive asymptotic expansions for weak gravitational fields 

We consider a gravitating system S made of ./V separated bodies labeled a, b or c (with a,b, c=\, 
2, AO, each of which being a spatially compact body made of some material medium, so that 
the velocity field v [Eq.(2.8)] and the field of proper rest-mass density p* are unambiguously 
defined. It is for such a given system, in particular for the solar system, that the celestial 
mechanics based on some theory of gravitation has to be built and physically tested. Obviously, 
an agreement between observations and predictions based on arbitrary approximations would 
prove neither the correctness of the theory nor that of the approximation method. One therefore 
has to introduce some definite method of approximation, which should have been mathematically 
tested, if not in the theory of gravitation, at least in other domains. We choose the method of 
asymptotic expansions, which has been quite extensively tested in many domains (see e.g. 
Kevorkian & Cole [28], Awrejcewicz et al. [29], and references therein). Strictly speaking, a 
consistent asymptotic expansion may be defined and tested only if we have a family (S e ) of 
gravitating systems at our disposal, depending on a parameter £ which should take values having 
as an accumulation point. Hence, we should associate a such (virtual) family to the given real 
system S, with that system corresponding to a given, small value £q of the parameter. Of course, £ 
should characterize the weakness of the gravitational field. As is usually done, we shall associate 
the weak-field limit with a situation in which NG is a good first approximation (indeed an 
extremely good one, in the solar system). We shall try to introduce this link without logical 
circularity. As in Section 4, we consider that any matter is in the form of a barotropic perfect 
fluid. However, we do not need that the state equation, say p* = F(p), is the same for each of the 
N bodies: we merely assume (for simplicity) that F is the same function inside any given body; 
since the bodies are separated, it is unnecessary to mark this by an index a, etc. 

5.1 The weak-gravitational-field limit in Newtonian gravity 

In NG, there are five independent fields: p, the Newtonian potential U^, and the "flat 
space-time" velocity, u = dx/dT. The five independent equations for those fields are: Poisson's 
equation (of course, there is no distinction in NG between the proper rest-mass density p* and 
the density of rest- mass in the given reference frame, and we shall note both simply p): 



At/ N = -4%Gp, 



(5.1) 



the continuity equation: 



dp 



+ div(pu) = 0, 



(5.2) 



dT 



and Euler's equation, 
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du 



du 

— + (gradu).u 



= pVt/ N -Vp. (5.3) 



Let these equations be satisfied for a "Newtonian gravitating system" and, for any £ > 0, set £ = 
e/Eo (where £o will be defined hereafter) and define new fields thus: 

f/ Ne (x,r)^2f/ N (x,<^), p e (x,T) = ^p(x,^T), u e (x,T) = %u(x,%T). (5.4) 
Note that U Ne< = t/ N , Peq= P and u g() = u. Moreover, in defining 

F e (p) = ^F^~ A p) [thus F e (p ) ee £2 Fi(e~ 4 p )], (5.5) 
one gets from (5.4) 2 : 

Pe(x,T) = % 2 p(x,%T). (5.6) 

It is immediate to verify that the fields p e , u e , U^ e and p e satisfy also Eqs. (5.1) to (5.3). Let us 
introduce an adimensional field- strength parameter as Misner et al. [4] 

£o = (t/Nmax/c Y\ U Nm . dX = max { U N (x, T ); (x, T)& R 4 }. (5.7) 

(L^Nmax is assumed finite; this will be the case, e.g., for an "insular mass distribution" quasi- 
perodic in time; if necessary, one may restrict the time domain to a compact interval.) From 
(5.4)i, it follows that [(f/N£)max/c 2 ] 1/2 = £ . Thus, with any Newtonian system, we may associate a 
one-parameter family of systems. All systems S e have just the same geometry, only the time scale 
increases as £ 1 (accordingly, the velocity decreases as £ ), the density and the potential decrease 

2 4 

as £ , and the pressure as £ , as £ tends towards zero. This is the weak-gravitational-field limit in 
NG. (Yet the family extends to strong Newtonian fields, for £ is arbitrary.) This result explains 
order-of-magnitude estimates partially justified by virial theorems (see e.g. Will [5], p. 89) and 
allows to use them in a definite asymptotic framework. 

5.2 A convenient choice of units 

In expanding the field equations (e.g., of the scalar theory), it will be essential to know the 
order of the various fields. It turns out that, in adopting a specific variation of the units with £, the 
situation becomes very simple. The change of units affects merely the time unit, which is 

—1 2 

multiplied by £ for the system S e , and the mass unit, which is multiplied by £ . Note that G is 
not affected by this change. But the time T is changed to T = £T, hence the velocity to u' = £ 1 u 
and, in particular, the velocity of light becomes c' = £ _1 c so that, in the new units, the small 
parameter may be taken to be £' = 1/c' (with £' = £ if c = 1 in the starting units). The mass density 
is changed to p' = £ 2 p, and the Newtonian potential becomes C7 N ' = £ 2 U^. Finally, the pressure 
is changed to p' = £ 4 p. Therefore, we have in the new units: 

C/ N ' e (x, T ) = £ " 2 f/ Ne (x, T ) = £ " 2 £ 2 Un (x, £T ) = £<f 2 £/ N (x, £ ~ l T' ) = C/ N ,(x, T ), (5.8) 

and in the same way u' e (x, T) = Ui(x, T), etc. Thus, in the new units, all fields are independent 
of £, and they have precisely the values they have (in the starting units) in the particular system 
Si, characterized by £ = 1. The spatial, as well as the temporal derivatives, will hence also stay 
independent of £ in the new units, hence in particular 
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(5.9) 



[As noted by Hinch [30], the usual notation (p = O(e^) implies only that q> is of order lower than 
or equal to k, thus we adopt his notation (p = ord(£^) when <p is of order k exactly.] Note, however, 
that this situation in which the fields p e , u e , e and p e are just the same (up to a change of 
units) independently of e, is special to NG. 

5.3 Definition of the weak-gravitational-field limit in the scalar theory 
As was said at the beginning of this Section, the proposed definition for a weak- 
gravitational-field limit will consist in defining a family (S e ) of systems, from the data of the 
given gravitational system of interest. Since, as was seen in Section 4, the natural boundary-value 
problem in the scalar theory is the initial-value problem, we have to define a family of initial 
conditions. Let us first state the two conditions that should be fulfilled by the family (S e ): 

(i) The curved metric y should tend towards the flat metric y° as £ — > 0. 

(ii) The fields p*, 27, p and u should be of the same order in £ as in the weak-field limit of NG, 
hence should be ord(e°) in the varying units of §5.2. 

In view of Eqs.(2.3)-(2.4), condition (i) amounts to say that the dimensionless scalar field/ 
=f e tends towards 1 as £ — > 0. Therefore, we are led to assume an expansion 

f e = 1 + (/)£ k + 0(£ k+l ). (5.10) 

We must determine the integer k. To this end, we shall enter this expansion into the field 
equation (2.6). We first get for the spatial part: 

Af £ =£ k A$ + 0(£ k+1 ), 

because the differentiation with respect to the spatial variables does commute with the expansion 
with respect to e, under certain uniformity assumptions for the expansions [ensuring that the 
spatial derivatives of the remainder term in (5.10) are still 0(£^ +1 )]. The two foregoing equations 
are valid in the fixed units as well as in the varying units, because the length unit is unchanged. 
To achieve the determination of k, we adopt the varying units, such that the small parameter £ is 
proportional to £' =l/c'. But therefore, we have to precise whether the time variable is T = £T, or 
c T = £T/£' T. This is exactly the departure point between the PNA (time variable T ) and the 
PMA (time variable c T' ), as we shall see. Accordingly, in a PNA (resp. in a PMA), we must 
assume that the derivative, with respect to T (resp. with respect to c'T' ), of remainder terms like 
0(£^ +1 ) in (5.10), is still 0(£* +1 ). Under this assumption, we find that the term involving time 
derivatives in Eq.(2.6) is 0(£* +2 ) [resp. ()(£*)]. Hence the left-hand side of (2.6) is ord(£*). On 
the right-hand side, a is given by Eq.(4.2). Hence, we get a E = ord(£ °) by condition (ii), 
therefore the right-hand side of (2.6) is ord(e 2 ). Thus one must have k = 2 in Eq.(5.10). We may 
set (in any units): 

f e =\- 2V e /c\ V e /c 2 = ord(e 2 ) (5.1 1) 

(the coefficient -2 is there for convenience; the ratio by c makes V homogeneous to the 
Newtonian potential). Eq. (5. 1 l)i defines V e , the important point is the estimate (5.11)2- 
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Next, we recall a recent argument [18] showing that the time variable in the varying units 
must be T ', i.e., a PNA must be chosen. The first point is that the asymptotic expansions will 
have first of all to be valid at fixed values of the chosen time and space variables. Since we 
choose to formulate the expansions in the varying units, it is then not equivalent to assume 
expansions with x' = T ' as the time variable (PN expansions) or with x° = c' T ' (PM 
expansions), because the ratio x' °/x° = 1/c' = £/c is proportional to the small parameter: it is easy 
to convince oneself that a PN expansion can not generally be rewritten as a PM one, and vice- 
versa [18]. Moreover, in the varying units, the time variable is 7" = eT where T is the "true" time 
(i.e. that measured in fixed units), hence c' T' = cT is (proportional to) the true time. But since the 
true orbital velocities in the system S e vary like £, the true orbital periods are like £ ~ . Hence it is 
7", not c' T oc 7, which remains nearly the same, as £ is varied, for one orbital period of a given 
body (a). In the PMA, the time variable would be the true time 7 = T'/£ and one thus should 
have at fixed x and 7, if p e is ord(e°) in the varying units [as imposed by condition (ii)]: 

p' e (x,7) = p (x,7) + O(£) 

(where index means zero-order coefficient in the expansion). But this means that the position 
of the bodies is approximately independent of £ at a fixed true time, in contradiction with the 
assumption that the orbital velocities vary like £. Finally, the reason why the choice of 7' (resp. 
of c' T ' ) as the time variable in the expansions may be called a PNA (resp. a PMA) is that the 
first choice leads to Poisson equations for the gravitational potentials, whereas the second one 
leads to wave equations [see Eqs. (6.6)-(6.8) below]. However, we have just seen that the second 
choice is not allowed in the "near zone" occupied by the bodies. 

Now we define the family of systems, i.e. the fields f e , p e , u e , and p* e , from the data of 
one gravitating system of perfect-fluid bodies in the scalar theory: we assume that the fields /, p, 
u, and p* = F(p), are known (for all x<= M and in some interval I = [0,Ti] for 7). We provisionally 
come back to the invariable units for clarity, and we set [cf. Eq. (5.11)] 

V=c\l-f)/2, (5.12) 

V max = max { V(x, 0); xeM}, e = (V m Jc 2 ) m , (5.13) 

£ = £/£q. We define / e , p e , u e [and p* £ = F^pg) with F e given by (5.5)], as the solution (assumed 
unique for all xe M and in some interval [0,r2(e)] for T) of the initial value problem 

, ? V(x,0) df £ 2^ 3 dV 

/ £ (x,0) = l-2<^ 2 ^-fA -^(x,0) = --^ — (x,0), (5.14) 

p e (x, 0) = £V(x, 0), u e (x, 0) = £u(x, 0), (5.15) 

for the field equations (2.6), (3.8) and (3.10). Equation (5.14)i means that, in Eq.(5.11)i, 

V e (x,0) = £ 2 V(x,0). (5.16) 

This is the simplest way to get the estimate (5.11)2. Also, Eqs. (5.15) are the simplest way to 
ensure that the fields p e and u e are ord(e 4 ) and ord(e) respectively, as in the weak-field limit of 
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NG. Moreover, u e = ord(e) means an £ dependence of the characteristic time; if one combines 
this with Eq.(5.11), one expects the asymptotic validity of Eq.(5.4)i for V : 

V e (x,r)~£ 2 V(x,£T)as£^0. (5.17) 

In uniform conditions, one may differentiate (5.17) with respect to T. This leads to define 
d T f e (x, 0) by Eq.(5.14) 2 . Thus, the initial data (5. 14)-(5. 15) is justified by (5.11), which follows 
from conditions (i) and (ii), and by the wish that the fields p e and u e behave as in the weak- field 
limit of NG, i.e. by condition (ii). In other words, the minimal conditions (i) and (ii) imposed to 
the weak-field limit lead essentially uniquely to the initial data (5. 14)-(5. 15). Conversely, one 
expects that the initial data (5. 14)-(5. 15) ensures that equation (5.17) holds and that 

p e (x, T)~$ 4 p(x, %T) and u e (x, T ) ~ % u(x, £T ) as £ — > 0, (5.18) 

whence also, the state equation being p* e = F^pg) with F £ given by (5.5), 

p* e {x,T)~i; 2 p*{x,l;T). (5.19) 

Thus, we have associated with the physically given system S a one-parameter family of similar 
systems, also governed by the equations of the scalar theory. This provides a framework in which 
one might test the convergence of the asymptotic expansions derived hereafter. Once again, the 
essential point is that, in deriving asymptotic expansions, the existence of a one-parameter family 
(S e ) is implicitly assumed: it seems better to have at least a (very) good candidate for such family. 
Moreover, the initial data (5.14)-(5.15) will allow us to state initial and boundary conditions for 
the expanded fields (Section 7). Note that Eqs.(5.10) to (5.19) hold true in the varying units, if 
one also expresses in these units the velocity of light and the fields V, p, u and p*, corresponding 
to the physically given gravitating system (which makes those fields depend on £ ...). 

Since Eqs.(5.17)-(5.19) may be assumed to be valid, all fields, e.g. p e ,u e , V e and p* e , are 
ord(£°) in the specific varying units, in which £ 1/c. This means that our condition (ii), as well 
as [entering (5.17) into (5.12)] our condition (i) are indeed verified. It also makes the derivation 
of asymptotic expansions a routine affair. There is no theoretical difficulty in considering 
expansions with an arbitrary number n of terms but, of course, the explicit expressions become 
very involved as n increases. Since, in the solar system, the second iteration, i.e., that 
immediately following the Newtonian approximation, leads merely to minute corrections, we 
shall stop at this second iteration. As a last preliminary point before expanding the fields in 
powers of £, one has to ask a question about the parity of the expansions. We proved that the first 
two terms in the expansion of the scalar field are even powers of £ 1/c [Eq.(5.11)]. But a look 
at the basic equations for a perfect fluid in the scalar theory, Eqs.(2.6), (3.8) and (3.10), shows 

2 2 

that it is only the square 1/c that plays a role (thus, the actual small parameter is 1/c ), hence 
more generally one may assume even powers in 1/c. In Ref. 18, some test of this assumption has 
been performed: it has been checked that the first-order expansion in 1/c gives the same results as 
the zero-order expansion in 1/c - namely, both expansions are equivalent to Newtonian gravity 
(NG). However, an important comment is in order: polynomials in the effective small parameter 
(1/c 2 ), i.e. Taylor expansions, are merely the most obvious possibility for asymptotic expansions, 
corresponding to the most regular behaviour. It may well be that these simple expansions work 
only until some finite order. Thus, in GR, there is a breakdown in the approximations based on 
even-powers expansions, either at the order 1/c 4 or at the order 1/c 6 , depending on the gauge 
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condition [8]. May be one should assume then new expansion functions, e.g. odd powers in 1/c - 
as is done for describing the first radiation effects in the standard PNA, where an 1/c 5 power is 
introduced [14]. 



6. Expansion of all fields vs. expansion of the gravitational field only 

In this Section, we shall assume that, in the varying units, the fields admit expansions with T 
(formerly T ' : we omit the prime henceforth) as the time variable; this may be called the "PN 
hypothesis". We begin with an ordinary asymptotic expansion, then we show that the standard 
PNA is a different kind of expansion. 

6.1 Ordinary asymptotic expansion: all fields expanded 

Since we assume only even powers in 1/c, the independent fields are expanded as 
follows: 

/= 1-2 U/c 2 -2A/c 4 + 0(£ 6 ), (6.1) 

p=p + p l /c 2 + O(£ 4 ), (6.2) 

u = u +ui/c 2 + O(e 4 ). (6.3) 

(Henceforth, we omit the index £ forf e , etc.) These equations are written in the (specific) varying 
units. The coefficients of the expansions are the fields U, A, p , pu Uo, Ui. By definition of an 
asymptotic expansion, these fields are independent of £ or 1/c (in those varying units), and they 
are also independent of the order of the expansion: in the next approximation, we add a term to 
any expansion, e.g. we add p2/c 4 to Eq.(6.2), but the fields of the preceding approximations, i.e. 
poand pi, are left unchanged. In invariable units, Eq.(6.1) holds true [although U becomes ord(e 2 ) 
and A becomes ord(e 4 )], whereas Eqs.(6.2) and (6.3) change to 

p = £ 4 [p +£ 2 p l /c 2 + O(£ 4 )], u = e[uo+e 2 u 1 /c 2 + 0(e 4 )], (6.4) 

where po, pi, etc., have the same numerical values as in the varying units (for the corresponding 
values of the respective time variables), hence do not depend on £. Having expanded the 
independent fields, we get by Taylor's formula an expansion of the other fields, e.g. 

p*=p* + p* l /c 2 + 0(£ 4 ), p* = F l (po), p*i=Fi(p ).pi (6.5) 

[the state equation has been assumed to be, in invariable units, p* e = F e (p e ) with F e given by 
(5.5); in the varying units, this becomes p* £ = Fi(p e )]. 

Let x° be the time variable used in the expansions. (Thus x° = T actually, and x° = cT for a 
PM expansion, even though we know from §5.3 that the latter is incompatible with the 
assumption that the true orbital velocities vary like £.) We get from (6.1): 

Af= -^-AU-^rAA+oU 6 ), (6.6) 
c c 
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-T^,o,o + -t( Ao,o + 2 U,o 2 + 4UU A0 ) + 0(e 6 ) . (6.7) 



Since x = T, the field equation (2.6) gives with (6.6) and (6.7): 



2_ 

~2 



-AU + - 



<9 2 [/ 



-AA + - 



<9r 2 



+ 



4tiG 



a + ^+o(e 4 ) 



(6.8) 



where on the r.h.s. we have replaced a [Eq.(4.2)] by its asymptotic expansion, which will be 
deduced a few lines later from the expansions (6.1)-(6.3). Thus, Oq and <5\ do not depend on £, 
hence we may identify the powers to get the separate and exact equations 

AU =-4-kGo , (6.9) 

AA =-47iGa 1 +^ 2 f; /dT 2 . (6.10) 

These two Poisson equations imply an instantaneous propagation of gravitation, as it is 
characterized at the order n = 2 by the fields U and A. It is easy to convince oneself that the same 
occurs at any approximation, i.e. at any order n of the expansion, under the "PN hypothesis" x Q = 
T. It is also easy to verify that the "PM hypothesis" x = cT would lead to wave equations with 
propagation at speed c - but the PM hypothesis is untenable in the near zone (§5.3 and Ref. 18). 

The spatial metric (2.3) is expanded, using (6.1), to 

^ = V^ (y + o(e 4 ), fl^fy-^/^ + ofe 4 ), (6.11) 



where 

i U i U i 

tt« - ' 3 . (6.12) 

3 u <k u, k 

We introduce the density of rest mass in the preferred frame and with respect to the Euclidean 
volume measure dV= Vp° dv 1 dv 2 dy 3 . This density (denoted by poo in previous works, while dV 
was denoted by dV° instead [20-21, 26]) is defined by [20]: 

p =dmo/dV=p*y u /p. (6.13) 

Its expansion is given by p = po + p/c 2 + 0(£ 4 ), with 

po = p* , Pi =p*i + Po(u 2 /2+ U). (6.13bis) 

Using this expansion we get by entering the expansions (6.2), (6.3) and (6. 1 l)i into Eqs. (4.2) 
and (3.9): 

y/b = Ob = Po = P*o, Oi = Pi + po(u 2 /2 + 77 + U), yfi = a i + Po- ( 6 - 14 ) 

The expansion of the energy equation (3.10) is then found to be 
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d T p Q + dj(p UQ J ) = 0, 



(6.15) 
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d i\ Po 

{ A J J " il v z ; 

= -podrU +d T po. (6.16) 

The local equations of motion (3.8) are expanded to 

d T (p UQ i ) + d j (p uju j )= PqUj-Pqj, (6.17) 

d T (p Q u x l +y/ 1 w i ) + o , 7 (Vi W + Po u i u o J + Po u o u i J ) + i T jk Po u o J ' u q + Po u o Jd r k ij= 

= -p hi + 2 kijpoj + iffi U i + po A, i - 2Up U, i , (6.18) 



where 



i^jk = k ij,k + ^fjfcj ~ kjkj , k {j = Uh l ij. 



(6.19) 



Now Eqs.(6.15) and (6.17) are just the Newtonian equations applied to the zero-order expanded 
fields of the scalar theory: Eq.(6.15) is the continuity equation (5.2); with Eq.(6.17), it implies 
Euler's equation (5.3), thus 



<M"o') + M o,7 u o i = u ,i- 



PoJ 
Po 



(6.20) 



They imply the validity of the Newtonian energy equation: 
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•+ p G d T U = 0.(6.21) 



Using the latter, the second-order energy equation, Eq.(6.16), simplifies to 

d T p l + d j (p 1 u j + p u l j ) = 0. 



(6.22) 



The latter shows that the rest-mass is conserved at the PNA of the scalar theory: together with 
Eq.(6.15), Eq. (6. 22) implies that 



^ + div(pu) = 0( £ 4 ). 



(6.23) 



But the rest-mass is not exactly conserved for a perfect fluid. Indeed a creation/destruction rate 
is predicted by the scalar theory, in a variable gravitational field under a fluid pressure. For a 
weak field, this rate is ord(e 4 ) and compatible with the experimental evidence on mass 
conservation [26]. We also note that the equations for the first PN corrections, i.e. the equations 
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for the order 1/c , Eqs.(6.10), (6.18) and (6.22), are all linear with respect to the "PN fields" A, 
pi, lj/i, Ui. This is typical of a perturbation method. 



6.2 Standard PNA: gravitational field alone expanded 

If we apply to the scalar theory the standard PNA, as exposed by Fock [1], Chandrasekhar 
[2], Weinberg [3], Misner et al. [4], or Will [5], then we expand the scalar field/, Eq.(6.1), but 
we leave p and u "unexpanded". Then Eqs.(6.6) and (6.7) still allow to write the l.h.s. of the field 
equation (2.6) as the l.h.s. of Eq.(6.8). If we wish, we even may express the matter source o, 
Eq.(4.2), as a kind of expansion: 

a= Ob + oVc 2 + 0(£ 4 ), a = p, o x = p(u 2 /2+n + u), (6.24) 

so that the whole of Eq.(6.8) holds true with this new definition of Oq and Oi. But since p, u and 
II are unexpanded and hence (of course) depend on e, the coefficients of this "expansion", Oq and 
<7i, now do depend on £ (this is called a "composite expansion"). As a consequence, we may only 
rewrite Eq.(6.8) as the approximate and unseparate equation 
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-AA + 
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J 




(6.25) 



Clearly, this unique equation for the approximate gravitational fields of the two first orders, U 
and A, is not enough - the more so, as the same occurs for the dynamical equations, of course. In 
order to get a solvable approximation scheme at the 1/c level, one must postulate at least one 
separate equation including 1/c 2 contributions, e.g. 



AU =-47iGp+0(e 4 ). (6.26) 



Fock [1, §68 sqq.] uses implicitly the equivalent of this minimum postulate for the PNA of GR. 
Other authors postulate two separate equations [2-5], which is actually equivalent. However, 
Eq.(6.26) is not a consequence of Eq.(6.25). In fact, Eq.(6.26) is not generally true if p is 
assumed to admit a usual expansion p = p + p/c 2 + 0(£ 4 ): (6.26) leads then, unless pi = 0, to 
the impossible equality 

Af/+47iGp =ord(£ 2 )+o(e 4 ). (6.27) 

Now if the independent fields/, p, u admit asymptotic expansions, Eqs.(6.1)-(6.3), then p does 
admit an expansion, given by Eq.(6.13bis). Further, at the initial time T = 0, we may easily 
calculate its coefficients (6.13bis) as functions of the initial data (5.15). In the varying units used 
to derive the expansions, p changes to p' = £~ 4 p , so that (5.15)i becomes 

p' e (x,0) = p data (x) = £ "V(x,0) (Ve), (6.28) 

where the notations of Section 5 have been adopted, i.e., p(x,0) is the pressure field in the 
physically given gravitating system (that corresponds to the value £o for e) at time T = 0, in the 
starting (fixed) units. Let us now come back to the notations of Section 6, thus again omitting the 
prime and the index £ for p in the l.h.s., and using the index for ordering the coefficients of the 
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expansions. It follows from (6.28) that the asymptotic expansion of the pressure p is simply, at 
the initial time: 

p(x,0) = p (x,0) + Pl (x,0)/c 2 + 0(e 4 ), p (x,0) = p data (x), Pl (x,0) = 0. (6.29) 

Using Eqs.(6.5) and (6.13bis), we then get 

p (x,0) = FxCpdataCx)), p, (x,0) = [p (u 2 /2 + C/)](x,0). (6.30) 

Equation (6.30)2 definitely proves that pi ^ at the initial time, hence also in some interval 
[0,Tb] for T, and so Eq.(6.26) is definitely incorrect in the frame of the usual method of 
asymptotic expansions. From (6.25), one may only derive 

AU =-47iGp+0(e 2 ), (6.31) 

instead of (6.26). Although Eq.(6.31) [together with Eq.(6.25) and the corresponding 
approximate equations of motion] is sufficient to calculate the post-Newtonian corrections, it is 
Eq.(6.26) that is required in order to sum the Newtonian contribution and those corrections so as 
to get a total post-Newtonian acceleration, correct up to and including the 0(e 2 ) terms. In other 
words: if one contents oneself with the asymptotically correct Eqs. (6.25) and (6.31), it does not 
help to add the post-Newtonian corrections to the "Newtonian" (zeroth-order) contribution, for 
the latter involves an error which is of the same magnitude as those corrections. According to the 
method of asymptotic expansions, the equations of the zeroth order, Eqs. (6. 9), (6.15) and (6.20), 
are exact equations (as well as the equations of any order), but they govern the zero-order 
expanded fields, not the exact fields. 

We conclude that the standard PNA does not pertain to the usual method of asymptotic 
expansions. The latter method leads instead to Eqs. (6. 9), (6.10), (6.15), (6.18), (6.20) and (6.22), 
which differ from the equations derivable from the standard PNA in a crucial point: in the former 
system of equations, one has two independent fields po and p\ to characterize the rest-mass 
density in the second approximation, as well as two fields uo and Ui for the velocity. We do not 
mean that the standard PNA is incorrect: it might find a justification in a different framework 
from the present framework, which is the usual method of asymptotic expansions applied to a 
one-parameter family of gravitating systems built from the physically given system. Similarly, 
the fact that this same framework does not allow to develop an approximation scheme involving 
wave equations for the gravitational potentials and valid in the near zone, does not mean that the 
standard PMA is incorrect when one applies it in the near zone. Finally, we emphasize that the 
present results are still far from allowing rigorous numerical estimates of the error done when 
one substitutes the expanded fields and equations to the exact ones. At the present stage, one just 
may expect that the error done with the first PNA, as envisaged here, is equal to a "moderate 
number" X [0.01 < X < 100 ?] times the next power of the small parameter, thus Ae 4 - for the 
value £o of the small parameter that corresponds to the physically given gravitating system, which 
is such that £q 2 ~ 10~ 6 as far as the solar system is concerned. 



7. Boundary conditions for the expanded fields 
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The boundary conditions for the expanded fields should be deduced from, or at least compatible 
with, the initial conditions (5.14)-(5.15) that define the family (S e ) of systems. In §6.2, we 
already found initial conditions for the expansions of the fields p and p, and these are Eqs. 
(6.29)2-3 and (6.30)i_2. Similarly, rewriting (5.15)2 in the varying units, we find an initial 
condition for the expansion of the velocity field u : 

u (x, r = 0) = u data (x) = eo" 1 u £q (x, T = 0), m(x, T = 0) = 0, (7.1) 

where u Eq is the velocity field in the physically given system Se (expressed in the fixed units), 
whereas uo and Ui are the coefficients in the expansion (6.2). 

Now let us determine the boundary conditions for the gravitational potentials U and A. 
We adopt the explicit notation of Section 5 and begin by rewriting (6.1) in that notation: 

f' e (x,T')= 1 -2U(x,T')/c' 2 -2A(x,T)/c A + 0(£ 6 ), (7.2) 

whence (under the relevant uniformity assumption for the expansion): 

d T 'f' £ (x,T' ) = - 2d T U(x,T' )/c a - 2d T A(x,T' )/ c A + 0(e 6 ). (7.3) 

On the other hand, rewriting (5.14) in the varying units, we get (recall that ^ = £/£q and £ = clc 1 ): 

2 V(x,0) df' 2 1 dV 

f ' e (x,0) = 1- 4^ -^=f (x,0) = — ^(x,0) . (7.4) 
c' z £ Z oT c' l £ dT 

Identifying the powers in 1/c', we obtain the initial conditions: 

U(x, T = 0) = V(x,T = 0)/£ 2 , d T U(x, T' = 0) = d T V(x,T = 0)/£ \ (7.5) 

A(x, T' = 0) = 0, d T A(x,T' = 0) = 0. (7.6) 

(By the way, we note that, had we chosen an exponent n ^ 3 for % in (5.14)2, the orders in (7.3) 
and in (7.4)2 would not have matched.) However, the equations for U and A, Eqs. (6.9) and (6.10) 
[with Ob and Ci given by (6.14)], are Poisson equations, hence elliptic, thus the natural boundary 
conditions for U and A are spatial-boundary conditions at infinity. No such condition is provided 
by (7.5) and (7.6), except at time T' = 0. We must hence provide independently the spatial- 
boundary conditions at infinity for U and A. Since we want to describe a spatially-compact 
gravitating system, we expect that the usual conditions for the Newtonian potential are relevant 
and thus postulate: 

(Vr) U(x, T) =0(l/r)andlVf/(x, r)l = 0(l/r 2 ) as r =1x1 ^00, (7.7) 

(Vr) A(x, T') =0(l/r)andlVA(x, r')l = 0(l/r 2 ) as r =1x1 ^00. (7.8) 

Since the expansion (7.2) is in particular valid for the physical system of interest, characterized 
by £ = £0, Eqs. (7.7) and (7.8) represent a constraint imposed on that system - more precisely, a 
constraint imposed on the field /e (x, T ). In fact this seems to be the second constraint we have to 
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impose on system S Eq {the first one is the requirement that £o [Eq.(5.13)] be finite and indeed 
small}: even the spatial compactness is not necessary to state the initial-value problem (5.14)- 
(5.15) that defines the one-parameter family of fields. And the existence of the expansions (6.1)- 
(6.3), which has been postulated in this work, is likely to be a mathematical consequence of the 
definition (5. 14)-(5. 15). Thus, the imposition of the boundary conditions (7.7) and (7.8) has an 
important physical meaning, for it is precisely at this point that we choose to describe a "quasi- 
Newtonian" situation. 

8. Conclusions 

i) The standard post-Newtonian approximations (PNA) cannot be regarded as asymptotic 
expansions in the usual sense (§6.2). The same is true for the standard post-Minkowskian 
approximations (PMA): the method of asymptotic expansions leads to a scheme differing from 
the standard PNA and PMA schemes in that, according to the former method, all fields 
(including the velocity and the mass density) must be expanded. The resulting ("asymptotic") 
PNA scheme is thus, in a sense, more complex than the standard PNA and PMA schemes, but it 
has an obvious status from the viewpoint of the approximation theory (although mathematical 
work would be needed, e.g., to prove that the obtained PN expansions are indeed asymptotic to 
the investigated families of solutions of the exact equations - as seems likely). Furthermore, it 
leads to local equations that are linear with respect to the fields of the approximation following 
the first ("Newtonian") approximation. 

ii) In the scalar theory at least (but in fact also in GR), due to the fact that the relevant asymptotic 
approximation is a post-Newtonian one with instantaneous propagation, the propagation of 
gravity at speed c has no effect, not only at the order 1/c, but even up to and including at least the 
order 1/c 3 [since we have obtained expansions with a remainder term 0(l/c 4 )]. This is in contrast 
to what happens in models that a priori assume some "aberration formula" for the gravitational 
force [31-32]. That the model used by Van Flandern [32] leads to enormous effects of 
propagation, is due to its first-order departure from NG. It turns out that things go differently in 
relativistic theories of gravitation. The cancellation of first-order retardation effects in a such 
theory was first obtained by Poincare [33] (see also Logunov [34]): he actually imposed this 
cancellation condition, as a constraint on his Lorentz-invariant modification of Newton's 
attraction law. 

iii) In the investigated scalar theory, there is actually no non-Newtonian effect at all below the 
order 1/c . Expansions in even powers of 1/c have been postulated here, because in appropriate 
units the actual small parameter that enters the equations is 1/c , and they have been found 
entirely consistent. However, we emphasize that, if one postulates instead first-order expansions 
in 1/c, one still finds that NG is recovered up to the order 1/c [18]. 

iv) Since several mass densities intervene in the "asymptotic" schemes proposed, one may expect 
that the (monopole) mass of a celestial body (a) will not be characterized, in those schemes, by 
only one number, but by several - e.g. by two numbers Mo" and M{\ at the first PNA. This may 
be checked by an adequate integration of the present local equations of motion in the volume of 
the celestial bodies, that provides global equations for the mass centers. This task has now been 
completed and the results (to be presented in a forthcoming paper) show that the additional 
masses M\ can, however, be expressed in terms of the Newtonian fields of the time T = {i.e., 
the initial time in the initial-value problem considered in Section 5). Now the masses are much 
like adjustable parameters in celestial mechanics, and one must generally expect that the 
adjustment is theory-dependent (simply because the equations of motion for the mass centers are 
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theory-dependent). This means that the Newtonian masses M N a , which are obtained by an 
adjustment based on the first (zero-order) approximation alone, are indeed correct only for 
Newtonian celestial mechanics: they coincide with the "true" masses of the first approximation, 
Mq, only up to an error which is a priori of the same order as the second-approximation 
corrections [21]. Therefore, it would be incorrect to a priori dismiss the investigated scalar 
theory on the ground that it predicts preferred-frame effects in celestial mechanics. But it remains 
to see if this theory actually provides a more accurate celestial mechanics than does the 
Newtonian theory. 

9. Bibliographical note 

Preliminary results of this work have been presented in Ref. 18. The first version of the present 
paper was submitted in February 1999 to another Journal, and rejected. It was pointed out to the 
author that the construction of a family of similar gravitating systems, like (S e ), had been 
previously studied, within GR, by Futamase & Schutz [35]. (This was not known to the author.) 
They also use the "weak-gravitational-field limit in Newtonian gravity", called by them the 
"Newtonian scaling theorem", and they also expand the matter fields as well as the gravitational 
field. With this in mind, the main new results of the present paper are the following ones: 

i) The scalar theory is considered, instead of GR. 

ii) In Section 5 here, the construction of a family of gravitating systems is deduced from the data 
of a general gravitating system with compact support, assumed to obey the gravitational 
equations. In contrast, Eqs. (3. 13)4-5 of Ref. 35 mean that the spatial metric is extremely particular 
at time T = 0, and does not seem appropriate to describe a general gravitating system. Moreover, 
we have shown that the demand of a Newtonian limit, as defined by conditions (i) and (ii) in 
§5.3, leads in a rather compelling way to our definition (5. 14)-(5. 15) of the family of systems. 
We have discussed the relation between the initial- value problem that defines the family, and the 
boundary conditions imposed on the expanded fields. 

iii) We have introduced variable units (depending on the small parameter e), in which the 
derivation of asymptotic expansions is straightforward and elementary. 

iv) Due to the simpler structure of the scalar theory, it has been possible to precisely compare the 
"standard" PNA scheme with the "asymptotic" PNA scheme and to conclude that the standard 
PNA cannot be regarded as an asymptotic expansion in the usual sense. Although Futamase & 
Schutz [35], as well as Rendall [8], emphasized that their ("asymptotic") PNA equations differed 
from those of the standard PNA in that they expand also the matter fields (in contrast to the 
standard PNA where only the gravitational field is expanded), they did not examine in detail the 
consequences as to the compatibility between the "asymptotic" and standard PNA of GR. 
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